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SECTION  I 


INTRODUCTION 

One  type  of  simulator  used  to  simulate  a vertically 
polarized  electromagnetic  pulse  (EMP)  produced  by  a nuclear 
detonation,  is  basically  a conical  antenna  over  a ground 
plane.  Near  the  base  of  the  bicone,  the  structure  consists 
of  a solid  surface,  constant  flare  angle  bicone  which 
gradually  transitions  into  a wire  biconical  structure. 

While  most  of  the  current  flow  is  radially  directed  along 
the  wire  portion  of  the  bicone,  circumferentially  distri- 
buted wires  are  used  to  add  to  the  capacitance  per  unit 

\ 

length  as  a wave  travels  along  the  bicone,  The  density 
of  the  circumferential  wires  and  the  wire  bicone  angle  are 
then  varied  along  the  cone  so  as  to  make  an  equivalent 

bicone  characteristic  impedance  at  any  point  which  roughly  I 

approximates  that  at  the  bicone  input.  A tapered  resistive 
loading  along  the  structure  is  used  to  damp  out  the  high 
frequency  currents  along  the  antenna  so  as  to  minimize 

{ 

diffraction  from  the  edge  of  the  bicone.  The  present 
structure  also  has  a flat  Copcap  which  is  also  loaded. 

The  circumferentially-directed  wires  are  not  orthogonal 
to  the  radially-directed  wires  and  consequently  most  of  the  j 

circumferential  symmetry  is  destroyed.  A moment  method 
solution  which  would  take  into  account  all  the  wires  and 


i 


the  solid  surface  feed  region  is  thus  unduly  complicated 
for  assessing  the  basic  performance  of  the  structure. 

In  view  of  the  intentional  design  of  the  structure,  however. 


it  seems  reasonable  to  approximate  the  structure  by  a bi- 
conical  surface  whose  flare  angle  is  the  same  as  that  of 
the  actual  bicone  feed  region,  In  the  following  sections, 
then,  a moment  method  solution  procedure  is  used  to  obtain 
the  static  charge  distributions,  capacitances,  and  effective 
heights  of  bicones  of  various  flare  angles,  both  with  and 
without  topcaps.  Derivation  of  the  integral  equation, 
formulation  of  the  numerical  solution  p rocedur e, and  the 
numerical  results  follow  in  the  succeeding  sections. 


1 


SECTION  II 


FORMULATION  OF  THE  INTEGRAL  EQUATION 


1 


An  integral  equation  for  the  charge  on  the  bicone 
is  obtained  by  requiring  that  the  potential  on  the  bicone, 
as  produced  by  the  charge,  equal  the  driving  potential 
between  the  cone  surface  and  the  ground  plane.  Referring 
to  Figures  1-3,  which  define  the  appropriate  geometrical 

quant  it ie s,  one  notes  that  the  ground  plane  may  be  replaced 
by  the  bicone  image  on  which  a surface  charge  of  the  same 
distribur  but  of  opposite  sign  exists.  The  resulting 

rhe  top  biconical  surface  in  terms  of  the  charge, 


poten ' 
is  t 


47TE, 


2tt 


P„  (r  ) 
c c 


II 

. 0 0 

2 7T  L sin  8 n 

*// 


cc 


r sin  8„dr  d<J>' 
c 0 c 


ccJ 


c t 


R j 
c t 1 


r dr  d<J>  ' 
t t 


0 ‘ 


(la) 


for  observation  points  on  the  cone  surface  and 


4’ie, 


2tt  L 


J 1 tc  tc- 


| r^  sin  8 _ d r . d <J>  ' 


L 0 0 

2 tt  L s in  8 


/ / 

0 0 


0 

' f 1 
0 (fj 


c c I R + R 

t t 1 1 


r dr  d t ’ 
t t 


= V 


0 

(lb) 


for  observation  points  on  the  topcap. 


In  (1),  the  subscripts  c refer  to  quantities  defined 
on  the  conical  surface  while  subscripts  t refer  to  the 


topcap.  In  double  subscripted  quantities,  the  first  subscript 
refers  to  the  surface  on  which  the  potential  is  observed 
while  the  second  subscript  denotes  the  origin  of  the 
source  contributing  to  the  potential.  The  plus  and  minus 
superscripts  refer  to  source  points  on  the  cone  and  its 
image,  respectively.  If  the  cone  has  no  topcap,  only  (la) 
with  p equal  to  zero  is  used.  Since  the  fields  are  $- 
symmetric,  the  observation  points  are  taken  along  the  x-z 
plane  where  p=0.  The  various  distance  quantities  in  terms 
of  cone  and  topcap  coordinates  are  all  of  the  form 


± / ± 

R = A - B cos  4';  r,s 
r s r s rs 


c , t 


whe  r e 


, ± 2 '2  . , 2 a 

A =r  +r  - 2 r r cos  9. 
cc  c c + c c 0 


B = 2 r r ' sin^  0. 
cc  c c 0 


+ ? ’ ? 2 2 2 
A = r“  + r “ - 2r  Loos  9 _ +L“  cos  9n 
ctct+c  0 0 


B = 2r  r sin  9_ 
c t c t 0 


±2*2  * 2 
A = r + r “ - 2r  L cos  9 + l “ cos  8. 

tc  t c + c 0 0 


B =2rr  sin  S. 
tc  t c 0 


1 1 


2 ' ? 2 2 2 

r“  + + L4  cos“0q(1  T 1 ) ^ 


tt 


2r  r 
t c 


(3) 


The  linear  charge  densities 
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may  alternatively  be  used.  Since  the  charge  density  is 
independent , the  cp  integrations  may  also  be  performed  in 
(1).  These  involve  integrals  of  the  form 
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is  the  complete  elliptic  integral  of  the  first  kind.  Hence, 
using  (5)  and  (4),  (1)  can  be  written  as 


SECTION  III 


NUMERICAL  SOLUTION  PROCEDURE 

The  method  of  moments  [ i ] is  used  to  obtain  a numerical 
solution  to  the  integral  equation  (7).  Since  at  the  bicone 
feed,  the  linear  charge  density  approaches  that  of  the  infinite 
static  bicone,  whereas  the  surface  charge  density  is  in- 
finite there,  it  is  appropriate  to  expand  the  linear  charge 
in  constant  pulse  functions  on  the  conical  surface.  On  the 
topcap,  however,  the  surface  charge  density  approaches  a 
constant  at  the  center  of  the  cap  while  the  linear  charge 
density  varies  linearly  there.  Hence,  on  the  topcap,  the 
unknown  is  taken  to  be  the  surface  charge  density,  which  is 
expanded  in  pulse  functions.  Figure  4 shows  the  pulse 
expansion  scheme  used.  Subdomains  of  widths 
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are  used  on  the  cone  and  topcap  surfaces,  respectively. 
The  subdomains  are  centered  at  the  points 
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on  the  cone  surface  and  at 
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on  the  topcap  surface. 


Endpoints  of  the  subdomains  are  at 
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If  we  define  a unit  pulse  function 
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then  the  charge  expansions  may  be  written  as 
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where  the  quantities  Q and  Q are  the  linear  and  surface 
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charge  densities  in  the  center  of  the  m t h and  nth  sub- 
domains  on  the  cone  and  topcap,  respectively.  Substituting 
(15)  into  (7),  using  (4b),  and  enforcing  the  equality  on 
both  sides  of  (7)  at  the  centers  of  the  subdomains  results 
in  the  matrix  equation 
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The  elements  of  the  "capacitance  matrix"  are 
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When  p “ m an J when  q=n,  the  integrands  of  'i>cc  and  . 

respectively,  are  logarithmically  singular; 


These  terms  are  thus  handled  numerically  bv  subtracting  out 
the  singularity  from  the  integrand  and  adding  its  analytically 


evaluated  integral  as  follows: 
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The  resulting  integrands  in  (21)  are  non-singular  and 

aay  be  numerically  integrated  using  standard  methods. 

Solution  of  the  (p+q)x(p+q)  matrix  equation  (16)  yields  the 

charge  coefficients  Q and  Q . Recall  that  the  former  is 

cm  tn 

a linear  charge  density  quantity  while  the  latter  is  a sur- 
face charge  density.  Equations  (A)  give  the  formulas  for 
conversion  between  surface  and  linear  charge  density  quantities. 
From  (9)  and  (10)  the  capacitance  to  ground  and  effective  height 


are  calculated  as 


SECTION  III 


NUMERICAL  RESULTS 


In  this  section,  numerical  results  are  presented 
for  the  linear  charge  distribution,  capacitance,  and  ef- 
fective height  of  a cone  over  a ground  plane  as  the  cone 
angle,  0^,  is  varied.  Results  are  given  for  configurations 
both  with  and  without  a topcap.  For  convenience  in  simu- 
lator design,  the  data  is  presented  assuming  the  permit- 
tivity of  free  space  rather  than  normalizing  the  data  to 
a medium  independent  form.  However,  for  other  applications, 
the  data  can  easily  be  scaled  to  apply  to  isotropic,  homo- 
geneous media  with  permittivities  different  from  that  of 
free  space. 

In  Figures  5,  6,  and  7,  the  linear  charge  density  is 

plotted  for  various  cone  angles  ranging  between  0^=  2.5° 
and  9q=  85°.  Note  that  as  the  cone  length  L tends  to  in- 
finity, the  linear  charge  density  at  any  point  on  the  cone 
should  tend  to  that  of  an  infinite  cone  over  a ground  plane. 
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where  V ^ is  the  voltage  between  the  cone  and  ground  plane  or 
half  the  voltage  across  the  terminals  of  a bicone  structure. 
Since,  for  the  static  problem,  only  r e 1 a t 1 v e dimensions  are 
important,  L-°°  is  equivalent  to  r -*■0 . since  from  either  point 


1 8 


arj-ed 


charged  cone 


of  view,  r /L-*-0.  Thus  the  numerically  computed  results 
c 

should  approach  (24)  as  r ->0  and  this  limit,  also  shown  in 

c 

Figures  5,  6 and  7,  provided  a convenient  check  on  the 
numerical  results.  Note  that  for  narrow  cone  angles  (Figure 
5),  there  is  very  little  charge  on  the  topcap,  as  one  would 
expect.  For  large  cone  angles  (Figure  7),  note  that  the 
charge  on  the  topcap  approximates  the  increase  in  the  charge 
on  the  cone  without  a topcap.  The  total  charge  in  the  latter 
case  includes,  of  course,  the  charge  on  the  top  surface  of  the 
cone.  For  moderate  cone  angles  (Figure  6),  a substantial 
portion  of  the  total  charge  resides  on  the  topcap.  With  no 
t o pc a p , howeve r , the  edge  condition  [2  ] requires  a more  sin- 
gular charge  at  the  edge  of  the  cone.  Furthermore,  the  com- 
puted charge  in  this  case  is  the  sum  of  the  charge  densities 
on  both  sides  of  the  cone  surface.  The  net  result  is  that 
the  total  charge,  with  or  without  the  topcap,  is  roughly  the 
same  for  all  cone  angles.  This  is  strikingly  evident  in 
Table  1 in  which  is  tabulated  the  capacitance  of  a cone  with 
and  without  a topcap  for  various  cone  angles  normalized  both 
to  the  slant  height  and  the  vertical  height  of  the  cone. 

From  the  tables,  it  appears  that  the  addition  of  a topcap 
increases  the  total  capacitance  only  by  about  3 Z for  moderate 
cone  angles.  The  computed  capacitances  agree  fairly  well 
with  the  rough  estimates  used  in  [ 3 ] . 


Table  L.  Computed  Capacitance  of  a Cone  over  a Conducting  Ground  Plane 


In  Figures  8-11,  the  capacitance  of  the  cone  is  plotted 


as  a function  of  the  cone  angle,  with  and  without  a topcap. 

In  Figures  8 and  10,  the  capacitances  are  normalized  to  the 
slant  he igh t , whe rea s in  Figures  9 and  11,  the  capacitances 
are  normalized  to  the  vertical  height  of  the  structure.  Also 
shown  are  the  corresponding  capacitances,  C^,  that  would  be 
computed  assuming  the  charge  distribution  to  be  that  of  an 
infinite  bicone,  Eq.  (24).  Using  (24),  C is  easily  found  to 
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£n  cot 


The  excess  or  "fringing"  capacitance  is  then  just  C - C^. 

As  a matter  of  interest  and  as  a check  on  the  reasonableness 
of  the  computed  capacitances,  the  percentage  of  fringing 
capacitance  is  plotted  in  Figures  12  and  13  as  a function  of 
Sq.  As  would  be  expected,  the  percentage  of  the  capacitance 
attributable  to  fringing  is  smallest  for  cone  angles  near  0° 
and  90°.  Near  9^  = 40°,  the  fringing  capacitance  adds  70%  of 
C to  the  total  capacitance. 

Table  2 lists  the  effective  height  of  a cone  structure 
with  and  without  a topcap.  Note  that  since  the  linear  charge 
density  on  an  infinite  cone  is  constant,  Eq . (24),  the  com- 

puted effective  height  neglecting  fringing  would  be  Lcos9q/2. 
Since  the  fringing  fields  increase  the  charge  density 
near  the  cone  edge,  this  figure  is  a lower  bound  on  the  ef- 
fective height  but  is  approached  as  9 ^ tends  to  0°  and  90°. 
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Table  2.  Computed  El  fee-live  Height  of  a Cone  over  a Conducting  Ground  Plane 


tig  urn  14.  Effective  lieigUl  of  a cone  with  a topcap. 


c on  o 


This  can  easily  be  seen  in  Figures  14  and  15  which  plot 
the  effective  heights  as  a function  of  cone  angle.  Note 
that  the  addition  of  a topcap  increases  the  effective 
height  by  about  7 ",  maximum.  This  is  also  in  fair  agreement 
with  the  rough  estimate  of  [ 3 ] . 
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SECTION  IV 


CONCLUSIONS 

An  integral  equation  has  been  derived  and  numerically 
solved  for  the  static  charge  on  a conical  antenna  over  a 
ground  plane,  both  with  and  without  a topcap.  Capacitance 
and  effective  height  data  show  that  there  is  almost  negli- 
gible increase  in  the  capaci'tance  when  a topcap  is  added  to 
the  conical  structure  and  that  there  is  but  a slight  in- 
crease in  the  effective  height. 

Future  studies  should  concentrate  on  refinements  of 
the  model,  which  should  include  a more  accurate  modeling 
of  the  wire  cage  structure  and  include  the  circumferential 
wires.  The  effect  of  tapering  the  bicone  angle  should  also 
be  analyzed  and  a more  accurate  model  of  the  feed  region 
should  be  employed.  This  improved  model  should  yield  more 
accurate  data  for  the  design  of  biconical  structures  for 


use  as  simulators. 
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